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I. 


INTRODUCTION 


This report presents the analysis techniques used in analysis of state 
vector observability obtained from SST data and describes software and numerical 
tests. 

All facets of the program are idealized in that the Earth is represented 
by a point mass, there are no transmission delays at the satellites, and all 
measurements are instantaneous. These implications mean that the sofware runs 
quickly, has small core requirements, and has easily specified input. It also 
means in general that errors are smaller than in a real world problem with modeling 
errors. Hence the program can be used easily to determine feasibility. 

The conclusions of the numerical work indicate that for short arcs - less 
than 0.5 km - there is not complete observability. For arcs longer than 1 km 
errors are on the order of 100 m, however modeling errors at this arc length 
may become very important. 
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II. PROGRAM DESCRIPTION 


This program is designed to analyze state observability with SST tracking 
through the use of a simple, spherical earth program. 

All facets of the problem are idealized, in that the Earth is represented 
by a point mass, there are no transmission delays at the satellites, and all 
measurements are instantaneous. 
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2.1 Input : Input consists of: 


Coordinates of the Station: 

= Station Longitude (deg) 

= Station Latitude (deg) 

Coordinates of the Relay Satellite: 
a = Orbit Semimajor Axis (km) 
e = Orbit Eccentricity 
i = Orbit Inclination (deg) 

^ = Longtiude of Ascending Mode (deg) 

a) = Argument of Perigee (deg) 

<p = Initial True Anomaly, (deg) 

Coordinates of Tracked Satellite: 

Same as Relay Satellite. 


Parameters of the Data Arc: 

A = Interval Between Readings (sec) 
N = Number of Readings 


Print Controls 


IP1 = 1 

Print 

IP2 = 1 

Print 

IP 3 = 1 

Print 

IP 4 = 1 

Print 

IPS = 1 

Print 

IP 6 = 1 

Print 

IP 7 = 1 

Print 

IP 8 = 1 

Print 

IP 9 = 1 

Print 

IP10 - 1 

Print 

IP11 = 1 

Print 

IP 12 = 1 

Print 

IP13-= 1 

Print 

IP14 = 1 

Print 

IP 15 = 1 

Print 


P sensitivity matrix, 
p sensitivity matrix 
p information matrix 
p information matrix 
p projection matrix 
p projection matrix 
p covariance matrix 
p covariance matrix 
p correlation matrix 
p correlation matrix 
combined projection matrix 
combined covariance matrix 
combined correlation matrix 
inverse of p information matrix 
inverse of p information matrix 
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Visibility Criteria: 

CTAN = Cutoff angle to relay satellite (deg) 

CTDS = Minimum distance between earth and intersatellite LOS (km) 
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2.2 Data : The following parameters are in block data: 


tt = 3.1415926536 

a) = 0. 72921159E-4 rad/sec Earth's rotation rate 

6 3 2 

y = 398601. km /sec Earth's gravitational constant 
SL, = 6378. km Earth radius 
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In the desired inertial system, 
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This is 
and velocity 


done for each satellite and we denote 


by (r^,v^) and the tracked satellite 


the relay satellite position 
by ( r 2 » v 2 > • 
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2.3 Initialization : After input, the positions and velocities of the station, 

and both satellites are computed in an inertial system coinciding with the Earth- 
Center-Fixed (ECF) system at time zero. The station initialization is detailed 
in the description of subroutine STATN. 

To initialize the satellites, we first define the magnitude of the radius 
vector, 

a(l-e 2 ) 

ss — ■ — . m , r . , • 

l+e cos <j) * 


then the radial speed, 

r = e sin <j> » P = a(l-e 2 ) ; 

and then tangential speed. 


r$ = n/pp /r . 


In a cartesian orbit plane coordinate system with x through the ascending 
node and 2 along the angular momentum vector, we have 

Xq = r cos (<t>+ixi) 
y Q = r sin (<j>+aj) 



Xq = r cos (<f>+aj) - r<J> sin (<}>+u>) 
y^ = r sin (<(>+w) + r$ cos (<j>+w) 



In a cartesian, equatorial coordinate system with & through the ascending 
node, we have 


y e = y 0 C ° S i 
z e = y 0 s±n i 

*e = *0 

y e = y o COS i 

z = y_ sin i. 
e ■'O 
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2.4 Main Loop : The main body of code is executed N times, once for each time 

point. The purpose of the main body of code is the generation of the N by 12 
matrices and defined by: 

.th 

i row o 3,^(0) , V;L (0) ,r 2 (0),v 2 (0) 


.th f , M(t i> 

i ro» ot s 2 s^fO.Vj^COJ.tjCOJ.^CO) 


This is implemented in several steps. First, using subroutine STATN, the 
current position and velocity (r^Ct^), of the station are computed. Then, 

using subroutine DSAMM, the current position and velocity of each satellite 
(r, (t.), v..(t.)), (r„(t ), v_(t.)) as well as the state transition matrices 

-LlllZiZl 


¥ 


1 


3r^(0) ,v 1 (0) 


3r 2 (t i ) ,v 2 (t ± ) 
2 3r 2 (0),v 2 (0) 


are calculated. 


Then, in subroutine HCOMP, the matrices 

3p(t ± ) 

H 11 3r 1 (t i ) ,v 1 (t i ) 

apCt^ 

H 12 ar 2 (t ± ) ,v 2 (t i > 

9p(t.) 

H 21 r 1 (t i ) ,v 1 (t i ) 


3p(t.) 

H 22 9r 2 (t i ) , v 2 ( t i ) 


til 

The i rows of and S 2 


are then computed in the main body of code by 
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th 


row of s 


r th 


1 " fH u V *u r.i 


row of s 


■2 - !H 2 , 


21 V b 22 r s ] 
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2.5 Output : Because this program will be used to explore observability, there 

are a large number of output processes. 

The sensitivity matrices and can be printed. 

The information matrices 


"i ' s l s i 


and 


W 2 - s t 2 r l 1 s 2 


where and are the measurement noises on range and range rate, respectively, 
can be printed. 


The projection matrices 


P 1 = S i S 1 


P 2 - S 2 S 2 


and 


Si a- 1 


S 2 “ 2 1 



* 


can be printed. These matrices are very important for determining the observable 
space since they are identity matrices on that space. 


f f -j* *|* 

The "covariance matrices" W^, and = (W^+V^) can be printed. "Covariance" 
is in quotes here because singularity of these matrices does not imply perfect 
information. 

f f f 

The "correlation" matrices obtained from W^, and can be printed. 
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PROGRAM SST ( li"Pt »T • OUT' J i.*T • TApfc5=l?\'Pl>T , T APc6=00TpUT ) 

.C •: ••-:.r~tSS"*a~S.' T^tfc' STATION LONGITUDE' ( !»*•* Ot.6R£Eb) 

; . C*** *tA>--iS IS THc STATION L A T T TUi.’t (IN DEGREES) 


V-V/tA'S IS . THE L-UW61 T.UUt OF T H t; RELAY SATELLITE (IN DEGREES) 

- : e<MHr" '>- r IS ThE s> ^r-f'A-JC'P AXIS OF THE TRACKED SATELLITE OK'tsIT 

C ( i".J K iL<>F«£TFRS) 

L * * TTS IS THE ST APT I fab ANGLE OF THE TRACKED SATELLITE 

-t) r.tsprE 6 si ~ "" 

T?ELT IS 'THE 'MEASUREMENT INTERVAL 
. . N IS TtrE. • F : .! i-'i d L A OF ME A S UR EH E N T S ( >'j ) 1U0) 

.... - . — CFtHMOhr/T hiV A'L7 p IrWfcvZW • « * P E 

COO lOL/GNEI/TT ♦ TA , TS « XL AO » C (bit) . TIME * ** GELT . PO ( 6 ) . 


* SI ( 1 00 * 12 ) •52(100,1?) *H 1 1 (1*6) ♦ H 1 ? ( 1 , b ) • H? i ( 1 , t ) * 6?2 ( 1 * o') , 

- # — THr* in 9 fry 9 Tf*P2 ( T. A') ,W1 (100*12) * V? ( 1 0 0.12) * 

*-.<5*ri ( 12,1 2 ) ? a FLAG ()2) 9 ATE.OP (12) *SW2 (12*1?) . 

3* ’? 1 ( 6 ) • Gfi'G A • Tr- 1 A T ( 6 * 6 ) ». S T ( t> . b ) *XN!i*TVEC(t>) , C-f T (6.6) .1 ( 6 » 6 ) , T T R ( A < 6 ) , 
' CO p m O N iv C ( 7 2 ) , 1 • 1 1 ; 

C O ‘-1 0 A, EHU • A 1 J ( 3 . 3 ) 9 oJJ ( 3 » 3 ) • C 1 J ( 3 . 3 ) .013(3*3) « T HE T . be. T AM 1 ♦ 

•■■■■:« GbtJPGT ♦G2 tG3 .'GAs Ob * t 3a • G7 • F O » E 1 • F2 « f 3 *"F A * Fb ? Ffe »+' 7 » O0A *PS12(fc>»6) 

..--- ... SR 1. • SH2 .3ET A 12 * SO 1 ♦ S02 

• r,i>i MENS I ON , CumB ( 2 0 0 , 12) .Cl (200*12) 

u I ME -ms 1 ON mPSG (12) 

DIMENSION PO0(6).P10(6),P20(6) 

... _ ...... 1 ( AyE-) - 

.. . '_rT ■' i R — P I / ]_ H 0 *.cT0 
r.i '* - iiEkii ih STfrflUN COORDINATES 

i 10 )-TSS-, X t^S* PF- - ' 

O* PtA-J m RELAY COORDINATES 

H £ A V ( 5 • 10 ) A 1 .El » K I I * riUMl « RL U Ml . PHI 1 

•-€**: •rEA'I); I N TRACKED. “SOTELO I.TE COORD I. NA T tS 

' ;.. . • WEffD ( b , 1 0 ) A2 . £ 2 * H 1 2,. 3GM2 * RLOF‘2 * PH 1 2 

PgnU LF ARC TIMING . 

- -pe.Au ( b , ib')1jtL t . -r--- 

C** PE A U PRINT CONTROLS 

H E A O ( b 9 2 0 ) 1 R 1 ♦ I P 2 » I P 3 . I P A » I p b * I P 6 , I P.7 » i P rt * I P 9 . I P 1 0 » I R 1 1 * I P 1 2 , I P 1 3 * 
. . . . ........ j / r .... j*. p.j.4 , j-p-T b . TP 1 fe . 1 pi 7 9 1 Pi b . IP1 S * I P2L 

20. FORMAT (201 1.) ; 

C#^- PLAu OiOISE STviLwS 

■ - - RF-4 0 (• b , I o ) S I (rP-.-S 'FG p u 

C-s:--::- PpAj IN VISIBILITY CrITEkI A 
PtA:j (5, 10)CTAN,CTf'S 
C * » . *K£-AO IN -APR I Ok I SIGMmS 

. z-Jk E.AU (5* 12.) (APS© ('I ) * 1 = 1 * 12) 

12 FORMAT ( t>E 1 0,4) 

— ■ — Rc'ivi} - — TN— APC7'" Eh't/PG I : ‘J T S ' 

READ (5, 12) (TAPC(I) ,1=1.12) 

*5R ITc (6*50) TSS* yLapS * FF 

I T E ( 6 » b b ) A 1 , 1 1 • P 1 1 . HOh 1 » R L U M 1 , P p 1 1 
\ii RITE (6*6 0 ) A P * £2 * H I 2 * b-UP 2 * RL ( .'M2 . Ph I 2 
i3H ITc. ( 6 * 6b ) GEL T * N 
• — -*R TTfc‘( 6,1 1 1 ) S J ( '' * 3 1 hS! ' ■ 

-:P I TE (6, 1 12) CTA'Y.CTuS 

.vP ITS (6*113) ( APS-.* (1) ,1=1*12) 


H -10 



P I T c { 6 » 1 1 4 ) ( T k n C ( T ) *1 = 1*1?) 

T 11 - ■ Fiv K?«rl t bn. * * SI fSH * S I tffrO */ 2 t 20 '.«) 

115 ^ (5X* *CT*iri*CT OS^/SE?!) • B ) 

C . 

1 13 -* ) *'/ 6 f 20 .b/5x**APS’<-> (7- 

1 1 4 F oP !v; mT ( 5X * *T akC ( 1 --6 ) '*/hfc: 20 • d/bX * *TaPC ( 7 — 
50 FOw.'tMT ( * 1 *,<(••/, -::-TSS • <L*pS » FF »/3fc20 . E ) 

■sSr E.OHm.T (Sx VFR;r» 1 1 * .« 1 1 « fcOMl ,RLO« 1 « PH 1 1 */bt 

671 ■ ?F.ri H H’A T ( 5 X • *;ft ? y JR ? . p J ? » R 0 ?> ? • R L 0 M 2 » P H I ? *• / 6 t 

p S: .■■ . FORMAT ( 5. X . T *.h!*/fc20 .6 *15) 

C ; TI-p.L l'7?. TT 0 f' ; S - • - : 

FF = 1 , -FF 
C T N = C T A W * U T 1 > F 

' - P-nT* (l.-E r*rl ) 

, KH I 1 — PJH 1 1 * u T 0 R 
«~P./ ( 1 »+€l*C0S (Foi l ) ) 

#d=£-l-#ST>v f-PH 1 1 ") *S'^pT“( 

RRtfD— SwPT ( 2 .')U*P ) /R 
P H p 0 '! = Phil + k L. 0 H 1 0 TO h 
'•-7 - C=CuS ( Pr-iPO'-T) ~ 

s=5 lliT (FHPOM) . 

Pi ( .1 ) =R*C 

■--- i -~ * ‘=>1(4 T — R O* C : “F? P H F)7Vs-~- - 

pi (?) =w*s 

‘-‘1(5) = H' 0 S + 0 P P 0 C 
•• c=tr 0 S tfti-i-n-iyfO.fr) 

' S=S 1-0 (PIT 4ft TOP ) 

:.- •; . ,.p£ o.(3 ) PPltf 2 .IBKC •••... , . , 

,.^,4....- y =p-i ts-) fr>“ T • - - 

Pl (?) =P 1 (?) *c 

Pl (5) =F 1 (5) *C 

*•: j 0=€0 5 ( PT’u*hl *y.Tfcp. ) . 

• 7 S= S I N < HO h. 1 * O T.O P )"'•• , ■ ■ 

PI 0 .( 1 ) =-P;l < 1 )'*C-P If?) FS 

; P‘i 0 ( 4- : ) =P'i : (-4 ; ) 1 ■•F5 ) *S 

P 1 0 ( 3 ) =P 1 ( 1 )*S+P1 (2)*C 
Pl 0 ( 5 ) =P 1 ( 4 ) -»S+P 1 ( b ) * C 

- pRI It: ( b * bb >PT0 

■ ’ P= A ( 1 ,r-E 2 R-r?) 

..pptsaph i 2*0TOk 

- : n=p/ (- ] •. +F 2*COS ( PH I ? ) ) •• - - 

R0=c2*SlN(PHI2) *$Q»T (Z° ! U/P) 
pp*-il.)=SoP T (7 o i J * P ) / P 

— p KRM-ji’4 — pH I g— + P i— O.!'! 2 0 T ' ) P 
C=C05 (PriPOP) ' ' 

. ' s=si^ (PHFioib;) ;; 

p:*(i ^ .... 

P2 (4) =PU*C-PPHO*S 

P? (?) =P*S 

P2 ( 5 ) =r U*S+PPH| HH, 

C— COS ( P I 2* pT(iP ) 

'■'S.= SlO..(Pl 2.t.0TUk) 

- pp ippp y—p? — 

r ?0 ( 6 ) =P2 ( 5 ) -"S 
°? (?) =P 2 (?) --C 


12) -/be 20 . t ) 

1?) */6t20»b/*l*) 

20 . 5 ) 

20.5) 
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P? < b) =P2 ( b ) C 

"_ J ~ ' ' " )~ 

'- " 3 = S I =') ( ? • 0 P c Vf 0 T 0 r ) 

- pro ( i ) =P2 ( l ) *c-p? < ? ) *s 

*•• w — ^'^> ifxirp = p 2 . ( 4 r * e - v ; r ( 5 T *5 " — 

P2'i (.2 ) =P 2 ( 1 ) *S + “? <? ) *c: 
p?0 ( s ) =P2 ( 4 ) «•? + ->? ( 5 ) 

- ; ) H ?u ' -■ --• — • - - ~ 

- T « •= T S S * i;) T p R ~ 

■„ h'r vfsCC -- 

Cfr'S ( >ft A-r)- 1 •• •• - 

pnn ( i ) =* t * c ,■ j s < t ^ r l 

PfiO K~< ) =h-£ *S nv ( T :> )"»€!_ 

POO' ( 3 ) = R £. * S 1 N (XL A E ) 

POO (•<*.) =-POPfi> ,4*pOO ( 2 ) 

- -Pfrfr ( n ) =OOOE'*ROO t j ) '* 

h 0 f > ( b ) — 0 » 

P 0 ( 5 ) — 0 o 

- 7 -- ; k 1 -><. H >*= PCK )- E 3 ) - • 

NftP -0 

.DO 431; 1 = 1 t 12 

.. .■- ---• — ( pto--()-*-y J >"-f j 4 31 

j\j a P — g A p + 1 
IJO 4 3? K = 1 i 1? 

- - r.- 31 IN H P *!<:) = 0:v-Q, 

.. 4 32 " v£? {'MAR j'Kj,sy 
. .=•: . ?- . so (map « l ) = si ob /Aprs < 1 ). 

■ — --SSBA'-P-vp » l)-=s l 3 R-f V 2 P S 3 ( 1 ) • - 

431 CORTIMjfc 

I 4 4:0=1 

-I‘l =RAP • — 

• ' • P = 0 Vv' ; • 

C "• v 3 - ' ' , . 

- I . 3 HtP ' l -’ - •• 

0 -- “• 

uO 10 0 0 1 = 1 . '.4 

-“-24 T~i !*h~ J^UfcLT + TARC ( I-ARC ) 

•J= J+1 - . . 

J. E '( X 104E « L£ # T-A R C . ( I A fi C + T ) ) G 0 TO 25 

: i-A'RG=T &R-G +2 - - - — 

IF ( T mRC ( I ApC ) ,i,r . 1 .r_30) bO T 0 2b 
J = 0 

- 00 - TO. 24 

• 2 %-,. ' CONTINUE - - 

. ' ipf-BC z “ 0OLE :::5TATR-( Pfi ? p 0 0 y OMSA • T I Mfe ) 

•-“-A'--—':-- T03-A-H-0-.- - - 3'-. :• - ' . - 

0 mL L I ! S APf'* ( T i Et » T G r ‘ AK » P 1 0 • P 1 0 ( 4 ) « P 1 * p 1 ( 4 ) * 0 » 2 ) 
0 u 23 j J = 1 » n 


lid 23\ KK — 1 ? b - 

--R PS II ( ) =PSI2 CJJtKK) 

23 COOT I POP 

-C Attr Afw.-( - t ffsf «-T O" AE 9 P2vi » P2u ( 4 ) , P3 , P? { 4 ) «()»?) 

C A L L VIS ( P 0 <i P 1 9 P ? « LT A i" 1 * C T U S » I St.fc ) 

J f ( i 3 F t . £ 0 p 0 ) C~. ij T f',- 1 0 0 (i 
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11=11+1 

CALL' nuu.-iP ( Ff- ) 

UQ 7-0 K=l*f> 

... f \ H — !\ + 6 

— WMtK : )=-&■, CO • 

SI ( i 1 • K6 ) =0 . cO 
S2(U,K) =0 .in 

:i2,(ii‘«f>fa) -0*r o - ■ 

UO;7o L=1 *6 - 

; ■■SiiIi-rK)=Sl (Il +. Hil(l.L) * 

|-, rf h ) -S 1 (1 r » rt-b ) + H 12 (1 • L )■ 

S 2 ( i 1 « K ) — S 2 ( J 1 i F ) + *"• 2 1 C 1 « L ) 
S?< LI »K6) =S2 ( 1 1 *«\s ) + F2? ( 1 *L ) 

f.o-- 

. IF H 1 . fcu. 10 ft ) Hfi r fi 26 
10.0 0 ;. v : COMflKut 

. '... Cfj’“ ;t T li'-'ifF ----- - ■■■ 

C « -* 

; O* c.vO -imIM LOOk 
! . . 

S • *1 = 11 

j;. If (iPl.NL.i ) ‘30 TO 401 

1 - ~ f-fcl- 6 ' v2 0 0 ) (• (S*l~Hr J ) i .0— 1 » 6 ) -» 1 

wPITi (6,225) 

: > ; S I TE ( 6 9 2 0 u ) ( (FI (I , J) ,J=7.12) * 

- 401 <iQivf fivbfe — ■ t - 

IF (IP^*sE . 1 ) SL- - TO : 4 n 2 . 

. •>i2 1T£ ('6,250 ) . 

• : " --—4- f Y- £ (• b-, - 2 (»•(•» ) ( (Si Of 4 J ) 'V J= I vb ) . I- 
■’< PITS (6,225) 

-21 Tc (6,200 ) ( (Si ( j. . j ) . j=7 , 1 ? ) * 

- • - / •*■£ i Y 25\<F): : ■ 

10 - h .. 

15 . FO'wM'A-T.ifSO _ 

' 2 ' 0 -O“ : 1 --4:: 

22b 2 0 2 T ( 5 ( / ) ) 

250 h OP vC| ( 1 0 (/ ) ) 

4 02 - .-iQQn T ir-Ht. - ... 

' i X: 0 .; i-, 1 , wi 

. _ '>:o iifi 4^1 v J? v 

~ ^- 5 -pf f « -J j =S .1 ( l « d )-/S J-fjH 

52 < i , J) =S2 ( I . J) /516 k U 

1 d.j)=si ( r » j) 

---.-■■<42. ( I r J ) =S2 ( 1 , d ) 

■ • • -f. '.-C’Odd ( J 9 d )—S 1 ( J. 9 J ) 

; . , t.^rf ,Ll.(i J) = C.U 0 -H ( 1- , J) 

- — „ j_. . . L 

12 = Nil + I - r-. ; A’P 
L 0 ( 1 2 « J ) = S 2 ( l , J ) 

... ... _ c i ( 1.2 9 d) = <3Uf,h < i 2 , J ) 

11-6 , .OOwTT-NUe. - 

-■ ' i •- (123 . Xc a 1 )-{?U TO 410 
. ■ . — — . — C <- L L' f' ?“r j-iL T ( r '>‘ 1 e —■ } • S i 9 1 9 1 (i‘0 f 1 2 

C^LL SPHi.fT ( S- 1 } 

IF ( I PI 4 .H . 1 ) (-.0 r< 4 1 0 


P S 1 1 ( L. 9 K ) 

<- P? 12 (L.K) 

* PS 1 1 ( L * F ) 

* P S I 2 ( L » F ) 


=1 » Ml ) 
i = 1 • m 1 ) 

= 'T » i' T ) 

1 = 1* p 1 ) 
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:4iU~- 
40 5 


4 0 6 

4 12 ■ ' 


416 

-41/- 

4 13: . 

413 

4 20 . 


"4=6 n 1 -- 
467 

421 


C c L L O' I N v' 2 ( > 1? 1 « ” 2 , > *. P L b 6 » 4 T F 4 P , ] 2 , 1 2 « 1 2 * 1 . t — 1 0 < 
. ~t Ji.LL '-‘ iSp rO^Tl 501) 

I P (1 P4 * Stir. . 1 ) ’dO TO 411 

Call S-Ful 1 ( n 1 • S2 • S2 » so l « 1 0 0 • 1 ? ) 

( S 1 • 1 ) - ■ 

If ( 1-15. fvfc . 1 ) tpv r<> 411 

CALL H I At V 2 ( S'*; 1 , ~> 'i 2 , hF L a b « A T F.-4P , IP <> 12* 12*1 • t— -1 0 < 

i - - • 

^OU^JlfHtF -' Vo. 

? sClstlt ':<3I-W.v2.trS 1, 3~ 1 • T'FL Av? , a TF MP »-100*wl* 12,1. E— 1 0 * 
2 -f^ ; l-f'*FL4-6 , A?t^P> 10 0 • 01 ,12,1. E- 10 * 
1 4 ( jV b . N t .1 ) '.Hi 1 < • a ()5 
C 4 l L SFULT ( * ‘ 1 « S 1 « 1 , 3 1 , 1 0 0 i 12) 

■ c a t s » i r ■ ■ 

cd«T I HU-t ••' •“ 

. ,.I.F.( 1*76'". WE,. 1 ) 00 TO 40 6 

f-jj ^<P. ’f’-P, S'i- l , 1 00, 12) ..... 

call spkl'T (SMi ) 

COHTl'LUE 

rr -eA-LL~~sHt.lt T~(-HT ? S TvSl ) S-Kl ,100,12) 

"■ IF ( IPT; HE.' 1) So TO 412 
C-6 L L S'PRNTiS W 1 ) 

■ - 

00 415 1=1,12 
60=SQRT ( Sv> 1 (1,1) ) 

* . TiE'('S'Tr?ETi (J ,”0.) 00 'TP' "4 1 5 

UO 417’ J= 1,12 

IF ( J .•(?£* I ) Oh TO 416 

--^vtl (-l- vc ,)^S41Ti ', JIFFO — 

60 TO 4 17 

SOI (0,1 ) -SO 1 (J, I ) /so 

• • .C"O r lTX ^OC-- •; - • - •■ • 

SO 1 ( l , r ) =i . ; V 

, Co/vfTTiiifr. 

-e-ALL -SPFFT ( F-i)'- ■' - ' - 

Call SOUL T ( N l , $2 , S2 , S4 1,11(0,12) 

IF ( iPB.Ot . 1 ) SO To 420 

Ca.li_ sppivr is-.i'i ) 

; i F ( 1 p 10 o fTE i 1 ) 60 TO 421 
/-. u f » v 465 . 1=1 , 12 : ... ... 

' T ( Srtof’t T ,-.F )- ) '■ - - - "■ - -■ • 

IF ( SQ.tQ. 0 . 0 ) GO 'TO; 465 
u 0 467 0=1,12 
- 1 Ft O «-6p. IT 6#* TO 4t 6 — 

.••3$yfc ( I , 0 ) —Sr 1 1 (.1,-0 ) "/ So 

. .CO,J{> r ..4b7 r L- 'i- ■ 

COOT 1NHJ& 

S v 1 ( I , I ) = i . 

-CO-.-rT INUL - 
(CALL SPRivT ( Sw 1 ) 

COLiJITaUL 

: -V2=\« T+Tj-l-OAP : - 

C ALL 0 1 iv 1/ 2 ( COLA , Sw 1 , AP t„/ j .b , A TL F-H , 2 00 , i*>2 , 12, 1 .5 — 1 
It- ( IP1 1 .Nt . 1 ) bC« TO 423 


OJt T) 

be. 1 ) 

Ur T ) 
u r. T ) 


0 , L-t T ) 
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C Li. h>f- ! !LT ( 2 » CbMo • C 1 * S ! '- ; i 9 20 0 912.) 
SPrt^T t$<4l ) 

. . 423 /'CA;tt -SKULT ■(«>?. CHMt . COr-sb * $L*q »200 9 1 2) 

.. - I :jr- ( 1 M 1 2 < Nt > 1 ) foU I P 427 

— e^!:t ■ SP P iM T (5 1 ) 

427 lF(IP13,^t.l)STwK 

00 44b 1=1 9 12 

- • - -••• S?7~S- : JPT ( Sw i t r-, I ) T~ 

-If ( 5ft *£.<?• 0 . 0 ) (3 0 TO; 4 PS 
• ■■ 01) 4.97 J=l*12 

~ ; t'f TO 4vs — - 

SO I (I,J)=Swl (1 « J) /o '.; 

<30 TO 407 

, 49b ■ s.w i ( J ♦ i ) =$•« t (\yvir/su 

'.. 497 COfiTTNUt - 

• S:2l ; (it 1 ) = 1» 

.4-^5 — eoTrmvrrf - - 


C 4LL SPPi\ T ( S*.‘ 1 ) 
STOP . 
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2.7 Subroutine HCOMP 


Purpose ; To compute the sensitivities of the current range and range rate read- 
ings with respect to the current twelve-dimensional state. 

Analysis : Let 


(r 0’ 

V 

= (x 0 . 

V 

V x o 5 

y o* 

V 

(r l* 

V 

= (x 1? 

y r 


*1- 

V 

( r 2> 

V 

= (x 2 , 

y 2 * 

z 2 9 X 2 9 

y 2 * 

V 


denote the inertial position and velocity (state) vectors of, respectively, the 
station, the relay satellite, and the tracked satellite at the current time. 


The data readings at the current time are range. 


P = l r ! - V + ' r 2 " r l' = p 10 + P 21 


and range rate 


(v-l - v ) / (r - r ) (v - v ) * (r 2 - rj 

p = + 


10 


21 


p 10 + P 21 *" 


The sensitivity matrices are 


ap _ = H 

3r 1 ,v 1 11 


x r K o 


X 2~ X 1 


10 


21 


y r y o 


y 2 _y l 


10 


21 



0 


’2 1 

, 0 , 0,0 ; 

P 21 J 
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0 , 0,0 


. 9p = H 
3 r 2’ V 2 12 


X 2~ X 1 y 2~ y l 9 Vfl 
P 21 * P 21 P 21 


P 10 P 21 P 10 P I0 P 21 P 21 


1 

0 

1 

iH 

•>> 

y 2 y l 

p 10 

y r y o 

+ " 21 

y 2~ y i 

P 10 

P 21 

P 10 

P 10 

P 21 

P 21 ’ 

z r z o 

• « 
Z 2" Z 1 

P 10 

z r z o 

+ * 21 

Z 2 _Z 1 

P 10 

P 21 

P 10 

p 10 

P 21 

P 21 


ap 


H 2i 


h lll ’ h 112 ’ h 113 ’ 


9P = H 
3r 2 ,V 2 22 


Vi.!2iY!i 

P 21 P 21 P 21 


y 2 y l P 21 y 2~ y l Z 2 Z 1 P 21 Z 2~ Z 1 
P 21 P 21 P 21 * P 21 P 21 P 21 

h 121 ’ h 122’ h 123 
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SUBROUTINE HCGfcP (FF).i - 

' ’ ' ' COMMON/ 1 N VAL /P I 0MET7MI.1 ,7VV RET""* — ^ ” ~ 

COMMON/ONE/T T , T a , TS • XLAM , R ( 6 » 6 ) , TIME ,B,Ot'.LT»PO (6) * 

* SI (100,12) ,S2(1 00,12) ,H11 (1*6) ,H12( 1,6) ,H21 (1,6) ,H22<1,6) , 

* TMP 1(1,6), T NiP 2 (1, 6), W 1(100', 12)., M2 (100, 12), 

* SW1 (12,12) ,AFLA3(12) ,ATEMP(12) »SW2(12,12) , 

'*■ PI (6) ,0M6A, THAT (6,6) ,ST ( r 6, 6 ) * XNU » T VtC ( 6 ) ,QT (6,6) ,T (6.6) ,TTR (6,6) , 

* P2 ( 6 ) ' - " " " ' " 

EQUIVALENCE (PO(l) ,X0) , (P0<2) ,Y0) . <PO(3> ,Z0) , 

* (P1(1),X1),(P1(2),Y1),(P1(3),Z1), 

* ( P 2 ( T ) » X 2 ) T ( P 2 tET » Y 2 ) » ( P 2~( 3 ) , Z 2 ) » 

* (P0(4) ,XD0) , (PO(5)-*YD0) (P0(6) ,700) , 

* (PI (4) ,XU1> ,.<P1 (5) , YD1) , (PI (6) , Z01) , 

* (P2 (4)TX02) , (P2 C5) iYOP) , (P2 (6) «Z02) 

X10=X1-X0 

X21=X2-X1 

Y10=Y1-Y0 

Y21=Y2-Y1 

Z10=Z1-Z0 

221*Z2-Z1 ’ - - 

XD1 0=XD1-XD0 
X02 I -X.02-X0 1 

Y01Q=YDI-YDO 

Y021 = YL'2-Y01 

Z010=ZD1-ZDO ’■ J. 1 

- ZD 21 -ZD 2-7 DT ' ------ • - ■ ■ - 

R01= SORT ( X 1 0**2 + Y 10**2 + Z 10**2) 

R12= SORT ( X21**2 + Y21**2 + Z2l**2) 

SH2 1- ( X 1 0*XD1 0 + Y 1 0*YU TO + Z10*ZD10) /R01 
SH22= ( X2 1*XD2 1 + Y21*YD21 + Z21*ZD21) /R12 
WRITE(6,20)R01,P12,SH21,SH22 

20 FORMAT (4E20. 8//) - • 

DO 10 J=l,6 
HI 1 ( 1 , j)=0.E0 
H12 ( 1 , J)=0.E0 
10 CONTINUE 

Hll (1 ,1)=FF*X10/R01 - X21/R12 

HI 1(1,2) =FF*Y 1 o/RO 1- — Y2 1 /R 12 

Hll (1 ,3)=FF*710/R01 - Z21/R12 
HI 2 (1,1)= X21/R12 

H 1 2 ( 1 , 2 ) = Y21/R12 — 

H12 ( 1 ,3) = Z21/P12 

H21 ( 1 ,1 )=FF*XD10/R01-XD21/R12-FF*SH21/R01*X10/P01+SH22/R12*X21/R12 
H 2 1 ( 1 , 2 ) = F F * Y D 1 0 / R 0 1 - Y D 2 1 / R T2 - F-F * S H 2 1 / R 0 1 * Y 1 0 / R 0 1 + S H 2 2 / R 1 2 * Y 2 1 / R 1 2 
H21 ( 1 ,3)=FF*ZD10/P01-ZD21/R12-FF*SH21/R01*Z10/R01+SH22/R12*721/R12 
H21 (1,4) =H 1 1 (1,1) 

H-21 (1,5) =H1 1 (1,2) — - - 

H21 (1,6) =H1 1 (1,3) 

H22 (1,1)= X021/R12 - SH22/R 1 2*X?1 /R1 2 

- H22 (1 ,?>= YD21/R12- S H2-2-/ RT 2 *Y 21 7 PT 2— - - 

H22 ( 1 * 3 ) = ZD21/R12 - SH22/R 1 2*Z2 1 /R 1 2 
H 22 (1,4)= H12 (1 , 1 ) 

H2 2(1,5)= H12(l,2) 

H22H ,6)= H 1 2 (1,3) . .... ' . 

End ’V- — - --- ■ ■ ■ 
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2.8 Subroutine STATN 


Purpose ; To compute the inertial position and velocity of the station at the 
current time. 


Analysis : Let 


(r 


0’ 


V = 


(x, 


0 * A 0 > y 0 > 


V 


denote the inertial position and velocity of the station at the current time. 
The Earth-center-fixed (ECF) coordinates of the station are input in terms of 
latitude X g (spherical earth) and longitude 0g, so the initial ECF criteria 
coordinates of the station are 


Rg cos 0 g cos X g 


Rg sin 0 g cos X g 


*E 


sin X 


S ’ 


where is the Earth radius. 


Because the inertial coordinate system coincides with the ECF system at t = 0, 
the inertial cartesian coordinates of the station at time zero are 


x o ■ r e cos e s cos x s 

y o " "e sin 9 s cos x s 

0 - V • , 

z o " h sin x s 


.0 0 

x 0 = - “e y 0 


.0 0 

y 0 = “e X 0 


.0 n 

z o = 0 
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where oo is the inertial rotation rate of the Earth, 
e 

The inertial cartesian coordinates of the station at the current time, t, 
are then given by 


0 


0 


x n = x_ cos w t - y n sm a) t 
0 0 e J 0 e 


0 . ' , 0 
y 0 = X 0 Sln U e 1 y 0 C ° S “e 


0 

Z 0 = z o 


x = - o> y 
0 e ^0 


V„ = 0 ) x„ 
y 0 e 0 


= 0 . 
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SUBROUTINE ST ATN ( PC » POO » OMG A * T I ME ) 
UI MENS I ON PQ ( 6 ) • POO (6 ) 

WT=0MGA*TIME 
C=COS (WT) 

S = S I N ( WT ) 

PQ(1)=POO ( 1 ) *C-POO (2) *S 
PO (2) =POO ( 1 ) *S + PO0 (2) *C 
PO ( 4 ) =— 0M6A*P0 ( 2 ) 

P0(5)=0MGA*P0(1) 

RETURN 

END 


11-21 



2. 9 Subroutine DSAMM 


Purpose : To propagate the six- dimensional state vector and to compute the six 

by six transition matrix. 

Analysis : The state vector is propagated by the closed form, nonlinear, two- 

body solution. The transition matrix is computed by the closed form representa- 
tion for perturbations about the two-body orbit. 
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Ill, NUMERICAL STUDIES 


In these studies the measurement noise was given a standard deviation of 3m 
(range) and 1 mm/sec (range rate). 

Rosman was taken at -82.88 longitude 

35.20 latitude 

Mojave was taken at -116.89 longitude 

35.33 latitude 

ATS-6 was taken at -94° longitude in Series I. 

SERIES I: 

In this series a geosynchronous relay satellite was used with a tracked 
satellite in circular, polar orbit at about 840 km altitude. Experiments were 
made with station location and with the relative positions of the satellites. 

Case 1. Equatorial station directly under relay satellite, both on x-axis. 
Tracked satellite starts at south pole (-z-axis) , with velocity along y-axis. 
Tracked satellite was tracked for one period (*’6100 sec). Observability rank = 9. 

In this configuration, there is complete negative correlation between the 
relay satellite z-component and the tracked satellite x-component. From the 
projection matrix, we find that the observable combination of z^ and x^ is 

.028’ 

_-. 166. 

Thus a perturbation 


1 

> 

N 

M 


0.166 

_ AX 2_ 


0.028 


is unobservable. This perturbation ratio corresponds very closely to 


42164 

7213 
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(where and a ^ are the orbit semi-major axes), indicating that a rotation of 
the orbit planes about the line of nodes (y-axis) is unobservable. This accounts 
for a rank decrement of one. 


There is also complete correlation between the relay satellite z-component 
and the tracked satellite y- and z- components. From the projection matrix we 
find that the observable space is spanned by 


.99999982 
-. A2620828E-3 
L -. A3925539E-6 


-. A3925539E-6 " 
-.0010306110 
. 9999989A 


From this it follows that any vector 


1 

> 

N* 

I-* 

1 


1 

CM 

< 

= 

23A6. 267763 

1 

CM 

•N 

< 

1 


2. A18092003 

■ — 


is unobservable. This seems very much as if y ^ is directly unobservable. However, 
the units here are km and km/sec which means that the components are of about 
equal importance. In fact this perturbation corresponds very closely to a rota- 
tion of the orbits about the x-axis, that is to a change of 


Az- 


1 


Ay 2 

= 

- 4i 2. 



3.0388 
7213. 
7.A338 


sin i 


There is also correlation between y^, x^, and x^ with 


r 

rM 

•>> 

< 

1 


. 7292115 7E-A 

. 17630615E-3 

Ax^ 

= 

1. 

-. 12855708E-7 

l 

> 

Nj 

l_ 


-. 12855708E-7 

.99999997 

_ 
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being observable and therefore 


1 

rH 

<1 

1 


1 

Ax^ 

= 

- . 7292E-4 

Ax 2 


. 1763E-3 

- 


. _ 


being unobservable. This perturbation corresponds very closely to 


a l " 


"42164 

V 1 

= 

-3.0388 

V 2 


-7.4338 



m _ 


sin a 


i.e. a rotation of both orbits about the z-axis. 


Thus we see that small rotations of the axis system are unobservable for 
this case. 

Observability is the same with range, range-rate, or both. 

The variances in this case cannot be determined because the information 
matrix is singular. However the "variances" appearing in the generalized inverse 
of the information matrix are lower bounds for the recovered variances. These 
gave standard deviations as shown in Table 1. One of the principal observations 
from this table is that combining different data types having different correla- 
tion structure can improve statistics far beyond the expectation on the basis of 
the increased number of observations. 

These variances are unrealistically low, of course, because of pseudo inver- 
sion; they do show, however, the variances that can be expected on the observable 
subspaces. To gain an idea of actual variances, a priori standard deviations of 
100 km were put on satellite positions with the results shown in Table 2. These 
are all larger of course. They show that, while only a three dimensional sub- 
space is unobservable mathematically, as a practical matter, very little is 
observable. The following variables are unobservable in any realistic sense: 
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Variable 

Standard Deviation 

Range 

Range-Rate 

Combined 

X 1 

.047 km 

4.069 km 

.002 km 

y l 

.000 km 

. 000 km 

.000 km 

Z 1 

. 013 km 

.027 km 

.001 km 


.056 m/s 

.168 m/s 

.002 m/s 

y l 

.126 m/s 

.523 m/s 

. 006 m/s 


.192 m/s 

.816 m/s 

.013 m/s 

x 2 

.080 km 

.157 km 

.004 km 

y 2 

.000 km 

.000 km 

. 000 km 

Z 2 

.086 km 

. 201 km 

.005 km 

x 2 

.079 m/s 

.473 m/s 

. 004 m/s 

y 2 

.065 m/s 

.107 m/s 

.004 m/s 

z 2 

.046 m/s 

.028 m/s 

.004 m/s 


Case 1. Standard Deviations Without A Priori 

Table 1 
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Variable 

Standard Deviation 


Range 

Range-Rate 

Combined 

X 1 

.048 km 

4.065 km 

.002 km 

y l 

100.000 km 

100.000 km 

100.000 km 

Z 1 

98.568 km 

98.568 km 

98.568 km 


7.292 m/s 

7.294 m/s 

7.292 m/s 

h 

.127 m/s 

.522 m/s 

.006 m/s 

i i 

42.627 m/s 

42.634 m/s 

42.627 m/s 

X 2 

16.862 km 

16. 863 km 

16.862 km 

y 2 

100.000 km 

100.000 km 

100.000 km 

Z 2 

. 086 km 

. 201 km 

.005 km 

X 2 

17.631 m/s 

17.637 m/s 

17.631 m/s 

y 2 

.065 m/s 

.107 m/s 

.004 m/s 

z 2 

103.061 m/s 

103.061 m/s 

103.061 m/s 


Case 1. Standard Deviation Using A Priori 

Table 2 
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Relay satellite 


downrange (y^) 
crossrange (z^) 
radial velocity (x^) 
crossrange velocity (z^) 

Tracked satellite downrange (y^) 

crossrange (x^) 
radial velocity (z^) 
crossrange velocity (x^) 

Naturally, much of the problem is caused by the very poor geometry of the 
station and satellites. 
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Case 2. This case is the same as Case 1 except for an eccentricity of 0.012 on 
the relay satellite. This gave a rank of 10 because the combination (y^, x^, x^) 
(rotation about the z-axis) was recoverable. However, the variance on y^ is 
very large and the errors project onto x^ and ^ 2 * 
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Variable 

Standard Deviation 

Range 

Range-Rate 

S Combined 

X 1 

0.105 

km 

4.168 

km 

0.005 

km 

y l 

6330. 

km 

3148. 

km 

3133. 

km 

. 

2 1 

0.030 

km 

0.026 

km 

0.001 

km 

X 1 

4556. 

m/s 

226. 

m/s 

225. 

m/s 

y l 

0.916 

m/s 

0.525 

m/s 

0.047 

m/s 

2 1 

1.522 

m / s 

0.829 

m/s 

0.071 

m/s 

X 2 

0.177 

km 

0.157 

km 

0.008 

km 

y 2 

0.618 

km 

0.065 

km 

0.000 

km 

Z 2 

0.173 

km 

0.204 

km 

0.010 

km 

x 2 

11030. 

m/s 

548. 

m/s 

546. 

m/s 

y 2 

0.111 

m / s 

0.108 

m/s 

0.007 

m/s 

*2 

0.069 

m/ s 

0.028 

m/s 

0.005 

m/ s 


Case 2. Standard Deviations 
Table 3 
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Case 3. Rosman station. Tracked satellite starts at south pole (-z-axis) , with 
velocity along y-axis. Tracked satellite was tracked for one period (-6100 sec). 
Observability rank = 11 and 12. 

Standard Deviations for this run are shown in Table 4. They indicate how 
precarious observability is for this configuration. The unobservable subspace 
using range rate is the same as that noted in Case 1 as being equivalent to a 
rotation about the z-axis, i.e. 


1 

> 

I- 1 

1 


1 

Ax^ 

= 

-. 7292E-4 

Ax 2 


-. 1763E-3 


In Table 5 are shown the standard deviations for the same configuration with 
a priori standard deviations of 100 km on satellite positions. 

We see that the same variables as in Case 1 are still unobservable in any 
practical sense: 


Relay satellite 

downrange 

(yi> 


crossrange 

(z-j^) 


radial velocity 

( V 


crossrange velocity 

(z^ 

Tracked satellite 

downrange 

(y 2 ) 


crossrange 

(x 2 ) 


radial velocity 

(z 2 ) 


crossrange velocity 

(x 2 ) 


Results tracking from Mojave were qualitatively the same. However, the change 
in geometry gave demonstrable benefits in the most unobservable subspace (y^, x^, 
x^) . In these components, which have the largest variances, the errors were halved 
while the remaining errors were essentially unchanged (Table 4.2). Note that this 
improvement is essentially the ratio of the aspect angle of Mojave (22°) compared 
with the aspect angle of Rosman (11°). 
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Variable 

Standard Deviation 


Range 

[ Range-Rate 

Combined 

X 1 

1435. 

km 

75. 

km 

74.619 km 

y l 

395982. 

km 

.000 

km 

21436.028 km 

Z 1 

79462. 

km 

4375. 

km 

4363.806 km 

X 1 

28959. 

m/s 

4.582 

m/s 

1567.691 m/s 

y l 

94. 

m/s 

5.631 

m/s 

5.616 m/s 

£ 1 

1671. 

m/s 

88.475 

m/s 

88.237 m/s 

x 2 

13597. 

km 

745.593 

km 

746.542 km 

y 2 

3859. 

km 

202.687 

km 

202.141 km 

Z 2 

11. 

km 

.737 

km 

.735 km 

X 2 

69820. 

m/s 

.878 

m/s 

3780.034 m/s 

y 2 

7. 

m/s 

.458 

m/s 

.457 m/s 

Z 2 

3973. 

m/s 

208.657 

m/s 

208.095 m/s 


Case 3. Standard Deviations Without A Priori 

Table 4.1 
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Variable 

Standard Deviation 

Range 

Range-Rate 

Combined 

x i 

1453. km 

76. 

km 

75.675 km 

y i 

200071. km 

.000 

km 

10843.268 km 

z i 

80448. km 

4435. 

km 

4422.330 km 

k i 

14503. m/s 

4.727 

m/s 

786.007 m/s 

h 

93. m/s 

5.593 

m/s 

5.577 m/s 

z i 

1738. m/s 

91.810 

m/s 

91.559 m/s 

X 2 

13767. km 

758.729 

km 

756.628 km 

y 2 

4012. km 

210.318 

km 

209 . 744 km 

Z 2 

11. km 

.730 

km 

.728 km 


35268. m/s 

.897 

m/s 

1910.983 m/s 

> 2 

7. m/s 

.454 

m/s 

.453 m/s 

* 

Z 2 

4130. m/s 

216.516 

m/s 

215.925 m/s 


Case 3. Standard Deviation Without A Priori (Mojave) 

Table 4.2 
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Variable 

Standard Deviation 


Range 

Range-Rate 

Combined 

X 1 

4.976 

km 

4.514 km 

2.683 km 

y l 

99.981 

km 

100.000 km 

98.518 km 

Z 1 

98.285 

km 

98.430 km 

55.288 km 


7.634 

m/s 

7.303 m/s 

7.203 m/s 

y l 

1. 313 

m/s 

.539 m/s 

.512 m/s 

4 i 

42.784 

m/s 

7.276 m/s 

6.820 m/s 

X 2 

16.548 

km 

16.839 km 

9.327 km 

y 2 

94.937 

km 

16.735 km 

15.306 km 

Z 2 

.537 

km 

. 203 km 

.182 km 

X 2 

17.515 

m/s 

17.637 m/s 

17.335 m/s 

y 2 

.293 

m/s 

.108 m/s 

.098 m/s 

*2 

97.849 

m/s 

17.237 m/s 

15.771 m/s 


Case 3. Standard Deviations With A Priori 

Table 5 


III-12 










Case 4. The same as Case 3 except that a longer tracking arc was used, tracking 
for about five orbits (~9 hr). Observability rank = 11 and 12. 

The unobservable space using range-rate is still the z-axis rotation, 

Ay 1 1 

Ax x = -. 7292E-4 

A x 2 -.1763E-3 

The standard deviations for this arc are shown in Table 6. Notice that the errors 
are beginning now to be reasonable, although still far greater than the measure- 
ment noise. These errors are small enough however so that a run with 100 km 
position a priori causes no significant change in the Range or combined tracking 
results. 
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Variable 

Standard Deviation 

Range 

Range-Rate 

Combined 

X 1 

.043 km 

.012 km 

. Oil km 

f 

y l 

5.582 km 

. 000 km 

1.414 km 

Z 1 

.759 km 

. 883 km 

. 177 km 

h 

.411 m/s 

.001 m/s 

.104 m/s 

y l 

.001 m/s 

. 000 m/ s 

.0002 m/s 

h 

.083 m/s 

.023 m/s 

.022 m/s 

X 2 

. 130. km 

. 033 km 

.030 km 

y 2 

. 194 km 

. 054 km 

.052 km 

Z 2 

.001 km 

.000 km 

.0002 km 

X 2 

.983 m/s 

.001 m/s 

.249 m/s 


.002 m/s 

.000 m/s 

.0002 m/s 

Z 2 

.200 m/s 

.056 m/s 

.053 m/s 


Case 4. Standard Deviations Without A Priori 

10 hr. Arc 


Table 6 
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Case 5: A comparison was made between the recovery possible when the tracked 

satellite orbit is in plan view and when it is "end-on". Because of visibility 
constraints, this comparison was made on a half-orbit arc. The end-on orbit 
started the tracked satellite over the south pole with velocity along the x-axis. 

Overall recovery was better for the plan view orbit in the sense that: The 

largest variances (y^ and y^) were in the end-on recovery; only x^, z^, y^ and 
x^ had smaller variances in the end-on recovery; and the smallest variances 2 .^ 
and y^) were in the plan recovery. Downrange velocity and radial position of the 
tracked satellite had the lowest variances in both cases. 
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Case 6: In this and the following case, an attempt was made to discover the 

effect of increased relay satellite a priori. The runs were made with the tracked 
orbit in plan view, Rosman station, a 30 second sampling interval over 2790 
records (94 data points), and a priori standard deviations on the relay of 10m 
and 1 mm/sec. 

The standard deviations were: 



a 

X 

a 

y 

a 

z 

r 

26 m 

393 m 

159 m 

r-r 

5 m 

96 m 

35 m 


4 m 

88 m 

31 m 


From this we see that range-rate is a very much better data type than range, at 
this level of relay uncertainty. When the standard deviations on the relay x, 
y, and z were increased individually to 100 m, the following results were ob- 
tained. 

a = 100 



a 

X 

a 

y 

a 

z 

r 

74 m 

1698 m 

469 m 

r-r 

6 m 

96 m 

35 m 

comb 

5 m 

96 m 

34 m 


o = 100 

y 



a 

X 

a 

y 

a 

z 

r 

26 m 

434 m 

160 m 

r-r 

9 m 

127 m 

76 m 

comb 

6 m 

97 m 

54 m 
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100 


cr = 

z 



a 

X 

a 

y 

a 

z 

r 

29 m 

408 m 

162 m 

r-r 

18 m 

96 m 

35 m 

comb 

15 m 

89 m 

32 m 


These results show that for this tracking geometry, range- rate is a better 
data type and is less sensitive to relay satellite position errors. 

When the standard deviations on the relay satellite were increased at one time 
to 100 m and 10 mm/sec, the following results were obtained. 



a 

X 

a 

y 

a 

z 

r 

76 m 

1701 m 

480 m 

r-r 

20 m 

145 m 

85 m 

comb 

19 m 

144 m 

84 m 
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Case 7: This is the same as Case 6 except that the tracked orbit is in end-on 

view. For the reference run, the standard deviations were: 


25 m 
6 m 
6 m 


689 m 
80 m 
79 m 


69 m 
6 m 
6 m 


Again, it appears that range rate is a far superior data type than range. 
When the standard deviations on relay x, y and z were increased individually to 
100 m, the following results were obtained. 


a = 100 
x 


a 

X 

a y 

°z 

108 m 

762 m 

118 m 

9 m 

83 m 

7 m 

6 m 

80 m 

6 m 


a = 100 m 


a 

°y 

a 

X 

z 

26 m 

847 m 

82 m 

22 m 

434 m 

24 m 

14 m 

294 m 

17 m 


26 m 
18 m 
15 m 


690 m 
80 m 
79 m 


70 m 
6 m 
6 m 
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These results continue to bear out the superiority of range-rate data; 
however, they do not indicate that range-rate is less sensitive to relay position 
errors in this configuration. 

When the standard deviations on the relay satellite were increased at one 
time to 100 m and 100 mm/sec, the following results were obtained. 



a 

X 

a 

y 

a 

z 

r 

112 m 

1064 m 

127 m 

r-r 

41 m 

708 m 

39 m 

comb 

38 m 

688 m 

37 m 


Basically cases 6 and 7 indicate that for short arcs in this high- inclina- 
tion satellite orbit, range- rate is a significantly better data type than range. 
Further, there appears to be an overall tendency for range rate to be somewhat 
less sensitive to relay errors than is range data. 


HI-19 



Case 8: In this case we examined the effect of inclination of the tracked satel- 

lite orbit on tracked satellite recovery. For this purpose we used the conditions 
of Case 7 (end-on orbit) and gradually decreased the inclination. Figure 1 shows 
the standard deviation of the position recovery. The z-recovery becomes increas- 
ingly poor as inclination decreases and when the orbit is equatorial, sensitivity 
to cross- track components is zero. 

Notice that this inertial parameterization (rather than HCL) appears to be 
a natural one in that uncertainties in x and y are virtually constant while z 
grows. If HCL were used, y and z would change from cross track and radial at 
i = 90° to radial and cross track at i = 0°. 
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300m 


200m 

Standard Deviation 
of Tracked Satellite 
Position Recovery 


100m 


Figure 1 


0 


15° 30° 45° 60° 75° 90° 

Tracked Satellite Inclination 
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SERIES II: This series had ATS-6 and GEOS-C in various configurations. 


ATS- 6 was defined by 


a = 42164.18909 km 
e = 0.000387768 
i = 0.875423° 

8 = 120.909439° (257.894711°) 
a) = 109.629318° 

<f> = 129.398243° 


GEOS-C was defined by: 


a = 7213.103 km 
e = 0.001313909 
i = 144.871022 

0. = 153.618899° (290.604171°) 
a) = 93.839036° 


<t> = 190.876765° 
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Case 1: Sampling interval 300 sec. Arc length 10.4 hours. This called for 125 

data points, but because of visibility constraints collected only 81. 

Observability Rank = 12. 

The standard deviations using Rosman station are shown in Table 7. When the 
number of data points was increased while keeping arc length fixed, the improve- 
ment, particularly in the least recoverable components, was only slightly better 
than would be expected on the basis of n/n . Hence the relative magnitudes of 
the numbers in Table 7 appear to be approximately correct for this arc, even with 
continuous tracking. 

The standard deviations using Mojave station appear in Table 8. Using com- 
bined data, there is an improvement of nearly 50% as the case in Series I would 
lead us to expect. This improvement, however, is not uniform between range and 
range rate. Presumably this is caused by the aspect of GEOS in this case. 

The correlation structure in all of the cases up to this point has been very 
poor, with the y^ - x^ correlation (downrange distance- radial speed) being about 
-0.99999 from Mojave and about -0.999998 from Rosman. Some of the following cases 
are directed toward trying to remove this high correlation. 
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Variable 

Standard Deviation 

Range 

Range Rate 

Combined 

X 1 

.038 km 

.027 km 

. 009 km 

yi 

3.580 km 

348.592 km 

. 967 km 

z i 

.418 km 

.455 km 

. 090 km 


.264 m/s 

25.409 m/s 

.071 m/s 

*i 

.002 m/s 

.006 m/s 

.000 m/s 

1 

.052 m/s 

.019 m/s 

.015 m/s 

x 2 

.127 m/s 

15.084 km 

.033 km 

y 2 

.355 km 

24.393 km 

.099 km 

Z 2 

. 050 km 

.031 km 

.014 km 

X 2 

.321 m/s 

34.270 m/s 

.087 m/s 


.488 m/s 

50.251 m/s 

.131 m/s 

* 2 

.069 m/s 

.074 m/s 

.014 m/s 


Case 1. Standard Deviations - Rosman 
TABLE 7 
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Variable 

Standard Deviation 

Range 

Range Rate 

Combined 

X 1 

.036 

km 

.026 km 

.010 km 

*1 

2.172 

km 

390.291 km 

.517 km 

Z 1 

.467 

km 

. 104 km 

.087 km 


.154 

m/ s 

28.449 m/s 

.037 m/s 

h 

.002 

m/s 

.007 m/s 

.000 m/s 

z i 

. 066 

m/s 

.036 m/s 

.017 m/s 

X 2 

.151 

km 

16.812 km 

.033 km 

y 2 

.035 

km 

27.357 km 

.006 km 

Z 2 

.064 

km 

.026 km 

.015 km 

X 2 

.215 

m/s 

36.110 m/s 

.051 m/s 

> 2 

.350 

m/s 

56.250 m/s 

.084 m/s 

Z 2 

.072 

m/s 

.041 m/s 

.016 m/s 


Case 1. Standard Deviations - Mojave 
TABLE 8 
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Case 2: Sampling interval 10 min. Arc length 24 hours. Because of visibility 

constraints, only 96 points were collected. The standard deviations are shown in 
Table 9. We can see that the variances are approaching acceptable levels. In 
addition, the correlations have decreased. The y^ - correlation while still 
high (-0.99986) has improved considerably from (-0.99999). 

The achieved variances noted in Table 9 are, of course, unduly optimistic 
because of the presence of significant modelling errors. 
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Variable 

Standard Deviation 

Range 

Range-Rate 

Combined 

X 1 

. 008 km 

.008 km 

.000 km 

*1 

. 288 km 

118.658 km 

.067 km 

Z 1 

.080 km 

.149 km 

.024 km 

*1 

. 021 m/ s 

8.649 m/s 

. 005 m/ s 

h 

.000 m/s 

.002 m/s 

.000 m/s 

' z i 

.005 m/s 

. 005 m / s 

.001 m/s 

X 2 

.024 km 

5.134 km 

.006 km 

y 2 

. 026 km 

8.305 km 

. 006 km 

Z 2 

.006 km 

. 009 km 

. 002 km 

x 2 

.029 m/s 

10.984 m/s 

.007 m/s 

*2 

.039 m/s 

17.105 m/s 

.009 m/s 


.014 m/s 

.026 m/s 

.004 m/s 


Case 2. Standard Deviations, 24 Hour Arc - Rosman 

TABLE 9 
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Case 3: In another effort to increase the accuracy of recovery, the eccentricity 

of the relay satellite (ATS-6) was increased to 0.1 and the run condition of 
Case 1 repeated. The results are shown in Table 10. Comparing these errors with 
those in Tables 7 and 8 it can be seen that the increase in eccentricity does 
indeed aid the recovery. However, the most significant reduction is in the down- 
range recovery from range-rate only. Aside from this, it appears simpler merely 
to use the Mojave station. 
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Variable 

Standard Deviation 

Range 

Range-Rate 

Combined 

X 1 

.051 km 

.013 km 

.012 km 

*1 

2.532 km 

4.730 km 

.747 km 

2 1 

. 451 km 

.989 km 

.079 km 

% 

.168 m/s 

.307 m/s 

.050 m/s 

h 

.012 m/s 

.025 m/s 

.004 m/s 

h 

.091 m/s 

.048 m/s 

.027 m/s 

x 2 

. 148 km 

.333 km 

. 037 km 

y 2 

. 369 km 

. 420 km 

.109 km 

z 2 

. 045 km 

.042 km 

.015 km 

* 2 

. 232 m/s 

. 453 m/s 

.067 m/s 

* 2 

.291 m/s 

. 619 m/s 

.086 m/s 

z 2 

.172 m/s 

.197 m/s 

.043 m/s 


Case 3. Standard Deviations 
Table 10 


III-29 









Case 4: A question that arises very quickly is the effect of knowledge of the 

relay satellite. In this case, a priori sigmas of 10 m and 1 mm/s were assumed 
for ATS and a 10.4 hour tracking arc used from Rosman. The standard deviations appear 
in Table 11 and we see that GEOS errors are down to the 2m, 2 mm/s level. 

Along with this, the correlation structure is greatly improved. Also note 
that range-rate appears to be a better data type than range. 

This result indicated that much shorter data arcs could be used, so an arc 
of 3.5 hours at a sample interval of 100 seconds was attempted. This gave ex- 
cellent results with good correlation and maximum standard deviation of 6 m and 
6 mm/s (range only). Again range rate was a better data type, (3 m and 2 mm/s). 

Pushing this still further, a one hour data arc was attempted with a 40 sec 
sampling interval. This gave reasonably good results, standard deviations were 
210 m and 224 mm/s (range only). Range rate was a better data type (74 m and 
58 mm/s). Results with combined data gave 66 m and 52 mm/s. 
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Variable 

Standard Deviation 

Range 

' 

Range-Rate 

Combined 

X 1 

.001 km 

.003 km 

.001 km 

*1 

.008 km 

. 008 km 

.008 km 

Z 1 

.010 km 

. 010 km 

.010 km 

X 1 

.001 m/s 

. 001 m/ s 

. 001 m/s 

y l 

.000 m/s 

. 000 m/ s 

.000 m/s 

Z 1 

.001 m/s 

.001 m/s 

.001 m/s 

X 2 

.003 km 

.002 km 

.002 km 

y 2 

.003 km 

.002 km 

. 002 km 

Z 2 

.002 km 

. 001 km 

. 001 km 

*2 

.002 m/s 

. 001 m / s 

.001 m/s 


.002 m/s 

. 001 m / s 

.001 m/s 

*2 

.003 m/s 

.002 m/s 

.002 m/s 


Case 4: Standard Deviations - 10 m ATS A Priori 

Table 11 
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Case 5: Several runs were made to explore the effect on GEOS recovery of larger 

ATS errors. These runs are 1 hour arcs and thus can be directly compared with 
Case 4, where results were: 


range 
r - r 
comb 


When cr was increased from 10 m to 1000 m, the results became: 
x 

r 

r - r 
comb 


When a was increased to 1000 m, the results were: 


r 

r - r 
comb 


When a of ATS was increased from 10 m to 1000 m, the results were: 
z 

r 

r - r 
comb 


Because of the complicated geometry it is difficult to explain the effects 
of the individual components. 

When the relay apriori was increased from 10 m and 1 mm/sec to 100 m and 
10 mm/sec at one time, the results were: 


182 

m 

215 m 

28 

m 

153 

m 

75 m 

73 

m 

27 

m 

75 m 

23 

m 


171 

m 

361 m 

53 

m 

10 

m 

238 m 

47 

m 

90 

m 

158 m 

31 

m 


172 

m 

218 m 

24 

m 

9 

m 

207 m 

41 

m 

9 

m 

69 m 

14 

m 


a 

a 

a 

X 

y 

z 

165 m 

210 m 

24 m 

9 m 

74 m 

15 m 

9 m 

66 m 

13 m 
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r 

174 

m 

690 m 

133 

m 

r - r 

21 

m 

229 m 

46 

m 

comb 

20 

m 

227 m 

45 

tn 


It is clear that range rate is a better measurement type than range, 
particularly when the ATS errors are small. In this configuration it does 
not appear that one can say that range-rate is less sensitive to ATS errors 
than is the range data, (see Series I, Cases 6 and 7). 


III-33 




An Example of Equivalence 


Lack of complete observability means that there is some combination of 
position and velocity components which cannot be determined from the data. This 
implies that the state vectors leading to a given measurement sequence are not 
unique. That is, the satellite state vectors can 'be perturbed and still generate 
the same measurement sequence. 

To verify this property of unobservability, we took the orbits of ATS-6 and 
GEOS-C from May 2, 1975, 23 hr. 30 min. to May 2, 24 hr, 0 min. and used the 
idealized analysis program to develop a direction in 12-space which was unobser- 
vable over this half-hour arc. This direction is given by the following vectors: 


A pos^ = [70239.03 

38673.52 

20260.21m ] 

A ve^ = [ -2.912502 

5.043442 

-.2959665m/ sec] 

A pos2 = [11233.49 

6944.440 

3138.174m ] 

A veL = [ 1.646808 

-2.278966 

0 ] 


When this increment of state was entered in the real-word simulation model 
(not the idealized model), changes in sensor readings occurred of about 0.2mm/sec 
in range rate and about 135m in range. These differences were constant over the 
arc. When this increment was increased by a factor of ten, the sensor discrepancies 
jumped to about 13mm/sec and 15000m, indicating that the changes are nonlinear 
effects. 
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Appendix A 


Lemma: (A'A)^ = A^A'^ = A^A^' 


Proof: We show that it satisfies the four Penrose Axioms. 


3) (A'A)^A’A is symmetric 


A + A ,+ A'A = A + (AA + )’A = A + AA + A 

= A + A 


which is symmetric. 


-j* 

4) A'A(A'A) is symmetric. To show this, we recall that 


B + BB' = B 


and 


B ,f = B t * 


Then 


A'AA + A ,t = A*A ,+ = (A + A) ' 


which is symmetric. 


1) A' A(A' A) + A' A = A'A 

= A^AA'A from 4) above 
= A’A 


2 ) 


(A’A) + A , A(A’A) + 


(A'A) + 

A^AA^A'^ from 3) above 
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Appendix B 


In this program, the primary method of calculating covariances, etc. is 
via the pseudo- inverse of the sensitivity matrix, rather than via the pseudo- 
inverse of the information matrix. Because this is not a typical procedure, we 
would like to outline the rationale for its use. 

First, let us agree that in most instances, the technique is impractical 
since it requires the retention and processing of a matrix having dimensions of 
the number of measurements by the number of adjusted states. In addition, the 
treatment of a priori information is somewhat more complicated. 

However the inversion process itself is numerically better conditioned. This 
fact can be illustrated in a number of ways. 

First, consider a sensitivity matrix, S, which is N by n with N>>n. The 
information matrix, W = S'S, is n by n. It seems very reasonable that independence 
of n N- vectors will be easier to detect than independence of n n- vectors. 

Second, note that for any symmetric matrix S, the eigenvalues of W = S'S 
are the squares of eigenvalues of S. While this does not hold for arbitrary 
square matrices, there is a tendency for the largest eigenvalue to more than 
square and the smallest topless than square, thus more than squaring the con- 
ditioning number ( = max |— [). 

j 

Thirdly, some examples will illustrate the effect. 

Example 1: Consider the nearly singular matrix 



The eigenvalues are 

X 1 = 1 + n/I+T = 2 
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If gaussian elimination is used on this matrix then the second step finds the 
matrix 



and the pivot element is -e. The angle between the column vectors of this matrix 
is approximately 

. -1 
a = sin 

(This criterion is of interest when using a Gram-Schmidt procedure for inversion, 
and is particularly attractive because it can be applied to nonsquare matrices.) 




If gaussian elimination is used on this matrix then the second step finds the 
matrix 



2 

£ 

and the pivot element is — . 
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This means that results which are obtained from single precision operations on 
S can be achieved using W only by forming W and operating on it in double precision. 

The angle between the column vectors of W is approximately 

. -1 e 2 

p = si n -g- , 

which leads to the same conclusions. 

Example 2: Consider the 2Nx2 sensitivity matrix 

” 1 1+e " 

1 1-e 

S = = 

1 1+e 

. i 1-e _ 

The angle between the column vectors is approximately 

. -1 

a = sxn e 




When the information matrix 


2N 

2N 

2N 

2N(l+e 2 ) 


is formed, the angle between the vectors is essentially squared and the inversion 
difficulty is increased. 
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